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ABSTRACT 

J.  A.  Tomlin  published  a  paper  [1]  on  meeting  required  multi- 
commodity  network  flows  at  minimum  cost.  He  formulated  this  problem 
in  both  node-arc  and  arc-chain  form.  The  node-arc  linear  program  was 
attacked  by  the  Dantzig-Wolf e  decomposition  principle  by  expressing 
the  derived  master  program  as  convex  combinations  of  the  extreme  points 
of  the  derived  subprograms.  In  this  note,  it  is  shown  that  this  prob¬ 
lem  is  really  a  special  case  of  the  problem  where  one  is  attempting  to 
meet  minimum  cost  multicommodity  flows  without  flow  requirements  on  the 
individual  commodities.  Tomlin's  algorithm  is  then  modified  to  solve 
this  more  general  problem.  When  this  is  done,  the  subprograms  are 
homogeneous  and  the  master  program  is  a  nonnegative  combination  of 
their  independent  solutions. 


• Ar 

Any  views  expressed  in  this  paper  are  those  of  the  author.  They 
should  not  be  interpreted  as  reflecting  the  views  of  The  RAND  Corporation 
or  the  official  opinion  or  policy  of  any  of  its  governmental  or  private 
research  sponsors.  Papers  are  reproduced  by  The  RAND  Corporation  as  a 
courtesy  to  members  of  its  staff. 
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INTRODUCTION 

Consider  a  network  [N, A]  with  n  nodes  and  m  directed  arcs  joining 
pairs  of  nodes.  Assign  (i,j),  the  arc  directed  from  i  to  j ,  a  capac¬ 
ity  b  ^  0  and  a  cost  c..  per  unit  of  flow.  The  c  's  are  such  that 

}  ij  J 

the  total  cost  on  any  directed  cycle  is  nonnegative.  Commodity  k, 
k  =  1,  ....  q,  is  identified  by  its  source  sk  and  its  sink  tR.  It  is 
required  to  find  a  multicommodity  flow  that  satisfies  the  capacity  con¬ 
straints  at  minimum  cost. 

Tomlin  [1]  published  a  paper  in  which  c_  ^  0  and  a  flow  rR  of 
commodity  k  was  required.  It  will  be  shown  later  that  this  problem 
is  a  special  case  of  the  one  treated  here.  Other  special  cases  treated 
by  this  more  general  formulation  are  maximizing  a  linear  combination 
of  the  k  distinct  commodity  flows  and  finding  an  efficient  routing 
when  the  value  of  a  routing  is  proportional  to  individual  commodity 
flows  but  this  value  can  be  offset  by  transportation  costs. 


. . .  a  ■ . 
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DECOMPOSITION 


Let  y.  .  be  the  flow  of  commodity  k  on  arc  (i,j)  and  v^  the  total 
flow  of  commodity  k.  Then  it  is  required  to  find  y„,  v^>  ar>b  m^n  2 


such  that 


(1) 


-E  =„£»?, 


ti.j) 


(2) 


(3) 


E 

k  =  l 


y  .  <  b .  .  a  1 1  ( i ,  j ) 

tj  ij 


“■E‘V’i>  + 

j=i 


- v,  for  i  =  s. 
k  k 


vk  for  i  ■  tk 


0  otherwise 


k  =  1 ,  •  •  •  i  q 


Letting  A^  be  the  node-arc  incidence  matrix  of  the  network;  d^  an  n 

vector  containing  -1  in  the  s^  position,  +  1  in  the  t^  position,  and 

0  elsewhere;  y,  the  vector  of  arc  flows  <  y. .  >  for  commodity  k;  b  the 
Jk  ij 

vector  of  arc  capacities;  c1  the  vector  <  c^j  >  °f  arc  costs;  and 
s  a  vector  of  slacks,  (1-3)  may  be  written  as 


(4) 


c'y.  +  c '  y„  +  ...  +  cy  =  Z  min 

1  2  q 


Iyl  +  Iy2  +  "  ’  +  Iyq  +  Is  =  b 


Vi +  dri 


*2y2  +  d2V2 


=  0 


=  0 


=  0 


Ay  +  d  v 
q-'q  q  q 
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,  4.  ,i  \>  =  Cl  mav  be  decomposed  into 

Note  that  solutions  to  +  dRvR  0  may  v 

inlt.  of  flow  along  paths  i»  %  “  Ck  and  al0,,B  CyClCS' 

;ince  all  cycles  have  nonnegative  cost  and  since  elimination  of  flow 
cycles  free,  a  feasible  solefion  cancel  yield  an  infeasible  solution, 

cycle  flows  need  not  appeal  in  an  optimal  solefion  and  can  be  eliminated 

,  w  1  be  the  set  of 

from  consideration.  Thus,  let  WR  -  •(  WR1>  *  ’  kNR  I 

solutions  corresponding  to  one  unit  of  flow  on  a  directed  path  from 

s  t0  t  .  W  is  then  a  set  of  points  that  span  the  set  of  solutions 

i/  k  ^  • 

.  Thpn  aDDlvinft  the  Dan tzig- 

.n  A  v  +  d  v  =0  which  contain  no  cycles.  Then,  applying 

\yk  k  k 

Uolfe  decomposition  principle  [2]  one  may  write  y,  as  a  nonnegatrve 
combination  of  tbe  elements  of  and  (4)  becomes 


N, 


N 


i 

VxXj  ( c  wii>  +  2^  qj  qj 

j"1  r 

vl  *q 

■c-'  /T  a,  4.  +V'\  (  I  W  .)  +  Is  =  b 

L>ij  V  +  •••  ^  u  « 


X.  .  2  0 

kj 

Ihe  number  of  variables  in  (5)  is  of  course  too  large  to  enumerate. 

„ever,  suppose  we  Have  a  basic  feasible  solution  and  let  ^  be  the 

.  i  1  ijni  i  o r-  for  the  row  containing  b  .  on  the 

orresponding  simplex  multiplier  tor  tne  r  lj 

,  ,  _  .  Tf  -  >0,  then  s. .  may  be  intro- 

ight- hand  side  and  n  -  <  ^  >•  nij  U 

need  into  the  basis  to  give  an  improved  basic  solution.  If  all  *  °. 

hen  any  solution  with 


(6) 


c  w.  . 
kj 


'  I  w  <0 
n  1  wkj  < 


will  yield  an  improvement.  If  no  w  satisfies  (6),  then  the  current 

kj 

solution  is  optimal.  If  arc  (i,j)  is  assigned  a  cost  of  c  -  TT^j> 

then  the  left-hand  side  of  (6)  is  equal  to  the  total  cost  of  the  path 

represented  by  w  .  Thus  the  search  for  a  solution  satisfying  (6)  may 
kj 

be  found  by  finding  the  shortest  path  between  s  and  t  for  k  =  1,  ...,  q. 

K  K 

Efficient  methods  for  finding  shortest  chains  are  given  in  [3-5].  Many 
of  them  require  that  the  network  contain  no  negative  cycles.  This  is 
assured  when  all  n. .  <  0  due  to  our  initial  restriction  of  the  c. ,s. 

The  phase  I  procedure  for  finding  a  starting  feasible  basis  to 


initiate  the  algorithm  is  accomplished  by  setting  Is  =  b. 
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pt  nw  RF91HREMENTS 

.  -♦».  ,,ac  reouired  to  minimize 

Tomlin  [1]  treats  the  problem  where  it  was  requir 

,  .  ,  rn  ,he  restriction  that  the  total  flow  of 

total  network  cost  subject  to  the  resen 

.  .  pClual  to  r  •  Furthermore,  it  was  assumed  that  all 

commodity  k  be  tqua  ^ 

.  1  of  the  problem  treated 

c  2  0.  This  can  be  treated  as  a  special  case 

J  .  _ nrH  firial  node  and 

mthis  paper.  For  each  k  on.  merely  ■««»..  “ 

,  it  to  s  (or  alternatively  to  it  from  t  ). 
directs  an  artificial  arc  from  it  to  sR  (or 

nf  r  and  a  very  large  neg- 
This  artificial  arc  Is  assigned  a  capaci  y  k 

r  slnce  none  of  these  artificial  arcs  belong  to  cycles,  our 
ative  cost.  bince  none 

initial  restriction  on  the  c_s  will  hoi 

■  m  differs  only  slightly  from  ours. 

The  node-arc  formulation  in  (D  differs  y 

Specifically,  (4)  is  replaced  by 


(7) 


,  ,  4 -  r ' v  =  Z  min 

:'y  +  c'y2  +  ...  +  c  yq 


iYl  +  Iy2 


+  ...  +  Iy  +  Is  =  b 


=  d. 


Vl 


V2 


A  yn 

qq 


=  d 


y  >0  al  i  k 
yk 


Here  d  is  the  vector  with  an  rR  in  the  sR  position,  a  -rR  in  the  tR 

^  1  .  all  other  variables  are  as  defined 

position,  and  aero  everywhere  else;  all  other 

before . 

This  problem  Is  also  treated  by  decomposition.  However,  the  sub¬ 
orns  A  y  -  d  ,  are  no  longer  homogeneous  and  consequently,  Instead 
programs,  AY  uq> 


iwuywj*  ■■gig' --it: 
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of  looking  for  a  nonnegative  combination  of  independent  solutions  to 
the  subprograms,  we  look  for  convex  combinations  of  their  extreme 


points.  Thus  letting  -  [w  ^ 
=  d^,  (5)  is  replaced  by 


w, ,  ]  be  the  extreme  points  of 

kN,  1 
k 


The  extreme  points  of  these  subprograms  are  also  paths  and  hence  the 

W,  defined  here  is  identical  to  that  defined  in  the  last  section.  In 
k 

form,  (8)  differs  from  (5)  only  in  that  one  must  add  a  convexity  con¬ 
straint  on  the  X  .  for  each  k.  However,  if  (5)  is  used,  one  must  also 
kj 

add  constraints  reflecting  capacities  on  the  added  artificial  arcs 
(i.e.  the  b  vector  in  (5)  is  larger  than  that  in  (8)  by  the  number  of 
commodities)  which  differ  only  slightly  from  the  convexity  constraints 


in  (8).  Thus  for  this  special  case  the  formulations  are  almost  identical. 


ARC-CHAIN  FORMULATION 


In  [1]  the  minimum  cost  multicommodity  problem  was  also  formulated 
as  a  linear  program  in  terms  of  its  arc  chain  incidence  matrix.  Spec¬ 
ifically,  rows  correspond  to  arcs  and  columns  to  chains.  The  program 
was  then  solved  by  a  simple  extension  of  a  method  due  to  Ford  and 
Fulkerson  [6]  which  has  since  been  extended  by  Wollmer  [7].  It  was 
shown  that  this  method  turned  out  to  be  equivalent  to  applying  de¬ 
composition  to  the  node-arc  program. 

In  this  paper,  we  will  not  go  into  the  details  of  this  method 
other  than  to  comment  that  if  the  more  general  problem  of  this  paper 
is  formulated  in  terms  of  its  arc-chain  incidence  matrix,  it  may  also 
be  solved  by  a  similar  extension  to  the  method  proposed  in  [6]  and 
that  this  method  is  also  equivalent  to  applying  decomposition  to  the 
node-arc  linear  program. 
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CONCLUDING  REMARKS 

1 1;  has  been  shown  that  the  problem  treated  in  [1]  is  a  special 
case  of  that  treated  here.  As  other  special  cases,  the  problem  treated 
here  includes  (i)  maximizing  a  linear  combination  of  the  individual 
commodity  flows  and  (ii)  finding  an  efficient  routing  that  takes  into 
consideration  the  value  of  both  the  individual  commodity  flows  and 
the  transportation  costs  involved.  The  former  of  these  is  accomplished 
by  attaching,  for  each  commodity,  and  artificial  node  and  an  artificial 
arc  directed  from  it  to  the  source.  The  artificial  arcs  are  given 
infinite  capacities  and  costs  whose  negatives  are  proportional  to  the 
linear  coefficients  of  the  commodity  flows  in  the  linear  function  that 
is  to  be  maximized.  For  the  latter  problem,  one  also  adds  these  same 
artificial  nodes  and  arcs,  the  new  arcs  having  infinite  capacity.  The 
cost  on  the  artificial  arc  for  commodity  k  is  the  negative  of  the  value  of 
a  unit  of  flow  of  commodity  k.  Thus  the  scope  of  problems  treated 
in  [1]  may  be  significantly  increased  by  relatively  small  changes  in 
the  problem  formulation  and  algorithms. 
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